XV.  On  the  Determination  of  the  Exterior  and  Diterior  Attractions  of 
Ellijisoids  of  Variable  Densities.  By  George  Green,  Esq., 
Cains  College. 


QRead  May  6,  1833.^ 

The  determination  of  the  attractions  of  ellipsoids,  even  on  the  hypo- 
thesis of  a uniform  density,  has,  on  account  of  the  utility  and  difficulty 
of  the  problem,  engaged  the  attention  of  the  greatest  mathematicians. 
Its  solution,  first  attempted  by  Newton,  has  been  improved  by  the  suc- 
cessive labours  of  Maclaurin,  d’Alembert,  Lagrange,  Legendre,  Laplace, 
and  Ivory.  Before  presenting  a new  solution  of  such  a problem,  it 
will  naturally  be  expected  that  I should  explain  in  some  degree  the 
nature  of  the  method  to  be  employed  for  that  end,  in  the  following 
paper;  and  this  explanation  will  be  the  more  requisite,  because,  from 
a fear  of  encroaching  too  much  upon  the  Society’s  time,  some  very 
comprehensive  analytical  theorems  have  been  in  the  first  instance  given 
in  all  their  generality. 

It  is  well  known,  that  when  the  attracted  point  is  situated  within 
the  ellipsoid,  the  solution  of  the  problem  is  comparatively  easy,  but 
that  from  a breach  of  the  law  of  continuity  in  the  values  of  the 
attractions  when  p passes  from  the  interior  of  the  ellipsoid  into  the 
exterior  space,  the  functions  by  which  these  attractions  are  given  in  the 
former  case  will  not  apply  to  the  latter.  As  however  this  violation 
of  the  lavl^  of  continuity  may  always  be  avoided  by  simply  adding  a 
positive  quantity,  for  instance,  to  that  under  the  radical  signs  in 
the  original  integrals,  it  seemed  probable  that  some  advantage  might 
thus  be  obtained,  and  the  attractions  in  both  cases,  deduced  from  one 
common  formula  which  would  only  require  the  auxiliary  variable  u to 
become  evanescent  in  the  final  result.  The  principal  advantage  how- 
ever which  arises  from  the  introduction  of  the  new  variable  u,  depends 
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on  the  property  which  a certain  function  V*  then  possesses  of  satisfy- 
ing a partial  differential  equation,  whenever  the  law  of  the  attraction 
is  inversely  as  any  power  n of  the  distance.  For  by  a proper  applica- 
tion of  this  equation  we  may  avoid  all  the  difficulty  usually  presented 
by  the  integrations,  and  at  the  same  time  find  the  required  attrac- 
tions when  the  density  p is  expressed  by  the  product  of  two  factors, 
one  of  which  is  a simple  algebraic  quantity,  and  the  remaining  one 
any  rational  and  entire  function  of  the  rectangular  co-ordinates  of  the 
element  to  which  p belongs. 

The  original  problem  being  thus  brought  completely  within  the  pale 
of  analysis,  is  no  longer  confined  as  it  were  to  the  three  dimensions  of 
space.  In  fact,  p may  represent  a function  of  any  number  a*,  of  in- 
dependent variables,  each  of  which  may  be  marked  with  an  accent,  in 
order  to  distinguish  this  first  system  from  another  system  of  s analo- 
gous and  unaccented  variables,  to  be  afterwards  noticed,  and  F may 
represent  the  value  of  a multiple  integral  of  s dimensions,  of  which  every 
element  is  expressed  by  a fraction  having  for  numerator  the  continued 
product  of  p into  the  elements  of  all  the  accented  variables,  and  for 
denominator  a quantity  containing,  the  whole  of  these,  with  the  un- 
accented ones  also  formed  exactly  on  the-  model  of  the  corresponding 
one  in  the  value  of  V belonging  to  the  original  problem.  Supposing 
now  the  auxiliary  variable  u is  introduced,  and  the  s integrations  are 
effected,  then  will  the  resulting  value  of  V he  a function  of  u and  of 
the  s unaccented  variable  to  be  determined.  But  after  the  introduction 


* This  function  in  its  original  form  is  given  by 
F = 


r p'dx  dy  dz 

{(y  - xf  + (y  - yf  -f  (z'  - z)2}T 

where  dx'dy  dz  represents  the  volume  of  any  element  of  the  attracting  body  of  which  p' 
is  the  density  and  x,  y,  z are  the  rectangular  co-ordinates ; x,  y,  z being  the  co-ordinates 
of  the  attracted  point  p.  But  when  we  introduce  the  auxiliary  variable  m which  is  to  be 
made  equal  to  zero  in  the  final  result, 

y. r p dx'  dy  dz' 

-I  {(x-  xf  4-  (y  - yf  -1-  (z'  - z)2  + m2}— 

both  integrals  being  supposed  to  extend  over  the  whole  volume  of  the  attracting  body. 
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of  ii,  the  function  V has  the  property  of  satisfying  a partial  differen- 
tial equation  of  the  second  order,  and  by  an  application  of  the  Cal- 
culus of  Variations  *it  will  be  proved  in  the  sequel  that  the  required 
value  of  V may  always  be  obtained  by  merely  satisfying  this  equation, 
and  certain  other  simple  conditions  when  p is  equal  to  the  product 
of  two 'factors,  one  of  which  may  be  any  rational  and  entire  function 
of  the  s accented  variables,  the  remaining  one  being  a simple  algebraic 
function  whose  form  continues  unchanged,  whatever  that  of  the  first 
factor  may  be. 

The  chief  objeqt  of  the  present  paper  is  to  resolve  the  problem 
in  the  more  extended  signification  which  we  have  endeavoured  to  ex- 
plain in  the  preceding  paragraph,  and,  as'  is  by  no  means  unusual,  the 
simplicity  of  the  conclusions  corresponds  with  the  generality  of  the 
method  employed  in  obtaining  them.  For  when  we  introduce  other 
variables  connected  with  the  original  ones  by  the  most  simple  rela- 
tions, the  rational  and  entire  factor  in  p still  remains  rational  and 
entire  of  the  same  degree,  and  may  under  its  altered  form  be  ex- 
panded in  a series  of  a finite  nuinber  of  similar  quantities,  to  each  of 
which  there  corresponds  a term  in  V,  expressed  by  the  product  of  two 
factors;  the  first  being  a rational,  and  entire  function  of  s of  the  new 
variables  entering  into  F",  and  the  second  a function  of  the  remaining 
new  variable  h,  whose  differential  coefficient  is  an  algebraic  quantity. 
Moreover  the  first  is  immediately  deducible  from  the  corresponding 
part  of  p without  calculation. 

The  solution  of  the  problem  in  its  extended  signification  being  thus 
completed,  no  difficulties  can  arise  in  applying  it  to  , particular  cases. 
We  have  therefore  on  the  present  occasion  given  two  applications 
only.  In  the  first,  which  relates  to  the  attractions  of  ellipsoids,  both 
the  interior  and  exterior  ones  are  comprised  in  a common  formula 
agreeably  to  a preceding  observation,  and  the  discontinuity  before 
noticed  falls  upon  one  of  the  independent  variables,  in  functions  of 
which  both  these  attractions  are  expressed;  this  variable  being  con- 
stantly equal  to  zero  so  long  as  the  attracted  point  ^9  remains  within 
the  ellipsoid,  but  becoming  equal  to  a determinate  function  of  the  co- 
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ordinates  of  p,  when  p is  situated  in  the  exterior  space.  Instead  too 
of  seeking  directly  the  value  of  V,  all  its  differentials  have  first  been 
deduced,  and  thence  the  value  of  V obtained  by  integration.  This 
slight  modification  has  been  given  to  our  method,  both  because  it 
renders  the  determination  of  V in  the  case  considered  more  easy,  and 
may  likewise  be  usefully  employed  in  the  more  general  one  before 
mentioned.  The  other  application  is  remarkable  both  on  account  of 
the  simplicity  of  the  results  to  which  it  leads,  and  of  their  analogy 
with  those  obtained  by  Laplace.  (Mdc.  C^l.  Liv.  iii.  Chap.  2.)  In  fact, 
it  would  be  easy  to  shew  that  these  last  are  only  particular  cases  of 
the  more  general  ones  contained  in  the  article  now  under  notice. 

The  general  solution  of  the  partial  differential  equation  of  the  second 
order,  deducible  from  the  seventh  and  three  following  articles  of  this 
paper,  and  in  which  the  principal  variable  V \s  a.  function  of  6 + 1 
independent  variables,  is  capable  of  being  applied  with  advantage  to 
various  interesting  physico-mathematical  enquiries.  Indeed  the  law  of 
the  distribution  of  heat  in  a body  of  ellipsoidal  figure,  and  that  of  the 
motion  of  a non-elastic  fluid  over  a solid  obstacle  of  similar  form, 
may  be  thence  almost  immediately  deduced ; but  the  length  of  our 
paper  entirely  precludes  any  thing  more  than  an  allusion  to  these  ap- 
plications on  the  present  occasion. 


1.  The  object  of  the  present  paper  will  be  to  exhibit  certain 
general  analytical  formulee,  from  which  may  be  deduced  as  a very 
particular  case  the  values  of  the  attractions  exerted  by  ellipsoids  upon 
any  exterior  or  interior  point,  supposing  their  densities  to  be  represented 
by  functions  of  great  generality. 


Let  us  therefore  begin  with  considering  p as  a function  of  the  .? 
independent  variables 

> ^2  > ^3  i 

and  let  us  afterwards  form  the  function 

dx^  dx-i  dxi dxl . p 


r=/- 


— (1) 


{ - xlj  + {x;  - Xif  + + (x,-  Xs'f  + 
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the  sign  / serving  to  indicate  s integrations  relative  to  the  variables 

.r/,  ar/,  xj,  and  similar  to  the  double  and  triple  ones  employed 

in  the  solution  of  geometrical  and  mechanical  problems.  Then  it  is 
easy  to  perceive  that  the  function  V will  satisfy  the  partial  differen- 
tial equation 

tPV  d^V  d^V  n — s dV 

“ = ^ + Sc?  + w + -iriu 

seeing  that  in  consequence  of  the  denominator  of  the  expression  (1), 
every  one  of  its  elements  satisfies  for  V to  the  equation  (2). 


To  give  an  example  of  the  manner  in  which  the  multiple  integral 
is  to  be  taken,  we  may  conceive  it  to  comprise  all  the  real  values 

both  positive  and  negative  of  the  variables  a;/,  xi, xl,  which  satisfy 

the  condition 


1 *JVti 

-4-  H — ■\- 


X.: 


+ 


(«) 


the  symbol  z,  as  is  the  case  also  in  what  follows,  not  excluding  equality. 


2.  In  order  to  avoid  the  difficulties  usually  attendant  on  integra- 
tions like  those  of  the  formula  (1),  it  will  here  be  convenient  to  notice 
two  or  three  very  simple  properties  of  the  function  V. 

In  the  first  place,  then,  it  is  clear  that  the  denominator  of  the 
formula  (1)  may  always  be  expanded  in  an  ascending  series  of  the 

entire  powers  of  the  increments  of  the  variables  x^,  x^, x,,  u,  and 

their  various  products  by  means  of  Taylor’s  Theorem,  unless  we  have 
simultaneously 

a;,  = x^,  X2  = X2, Xs  = Xs  and  u = 0; 

and  therefore  V may  always  be  expanded  in  a series  of  like  form, 
unless  the  ^ + 1 equations  immediately  preceding  are  all  satisfied  for 
one  at  least  of  the  elements  of  F'.  It  is  thus  evident  that  the  func- 
tion F'  possesses  the  property  in  question,  except  only  when  the  two 
conditions 
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^ J-  _ 4-  _ _i_ 

d.^  di 


H — 7^  /:  1 and  7i  = 0 
d 


(3) 


are  satisfied  simultaneously,  considering  as  we  shall  in  what  follows 
the  limits  of  the  multiple  integral  (1)  to  be  determined  by-  the  con- 
dition (a)*. 


In  like  manner  it  is  clear  that  when 


X,  Xo 
^ + ^ + 


a 


+ 


(4), 


the  expansion  of  V in  powers  of  u will  contain  none  but  the  even 
powers  of  this  variable. 


Again,  it  is  quite  evident  from  the  form  of  the  function  V that 
when  any  one  of  the  « -i-  1 independent  variables  therein  contained  be- 
comes infinite,  this  function  will  vanish  of  itself. 


3.  The  three  foregoing  properties  of  V combined  with  the  equa- 
tion (2)  will  furnish  some  useful  results.  In  fact,  let  us  consider  the 
quantity 


j'dx^  dx^...dxs du 


2 

+ 


where  the  multiple  integral  comprises  all  the  real  values  whether  posi- 
tive or  negative  of  x^,  x^, a;,,  with  all  the  real  and  positive  values 

of  7t  which  satisfy  the  condition 


/>•  2 /yi  2 

1 *V2 

-T  + — + 

a: 


(6) 


• The  necessity  of  this  first  property  does  not  explicitly  appear  in  what  follows,  but 
it  must  be  understood  in  order  to  place  the  application  of  the  method  of  integration  by 
parts,  in  Nos.  3,  4,  and  5,  beyond  the  reach  of  objection.  In  fact,  when  V possesses  this 
property,  the  theorems  demonstrated  in  these  Nos.  are  certainly  correct:  but  they  are  not 
necessarily  so  for  every  form  of  the  function  V,  as  will  be  evident  from  what  has  been 
shewn  in  the  third  article  of  my  Essay  on  the  Application  of  Mathematical  Analysis  to 
the  Theories  of  Electricity  and  Magnetism. 
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«!,  a,  and  h being  positive  constant  quantities;  and  such  that 

we  may  have  generally 

a,  > dr. 


In  this  case  the  multiple  integral  (5)  will  have  two  extreme  limits, 
viz.  one  in  which  the  conditions 


X, 


X^‘’ 


s H 1 + ' 

ar 


X,‘ 


u 


+ — ^ + ji—  1 M = a positive  quantity. 

G/g  fl/ 

are  satisfied;  and  another  defined  by 


•(7) 


^2  2 )|»  2 

+ -^+ + / 1 and  u = 0. 


a./ 


a, 


Moreover,  for  greater  distinctness,  we  shall  mark  the  quantities  be- 
longing to  the  former  with  two  accents,  and  those  belonging  to  the 
latter  with  one  only. 


Let  us  now  suppose  that  is  completely  given,  and  likewise  V/ 

or  that  portion  of  V in  which  the  condition  (3)  is  satisfied ; then  if 
we  regard  F"/  or  the  rest  of  F''  as  quite  arbitrary,  and  afterwards  endea- 
vour to  make  the  quantity  (5)  a minimum,  we  shall  get  in  the  usual 
way,  by  applying  the  Calculus  of  Variations, 

0=  — fdx^dxo dxsdu  w”"* ^ -h  ^7-^  + - — 

[ dx^  did  u dll  J 

dV 

— fdx■^  dx^ dxs  ^ V2  (8) 

seeing  that  ^ V"  — 0 and  S = 0,  because  the  quantities  V"  and 
are  supposed  given. 


The  first  line  of  the  expression  immediately  preceding  gives  generally 

...(2') 


^ ^ ^.d^V  d^V  n-sdV 

0 = + 


dxr  ‘ did  ' u dll 
which  is  identical  with  the  equation  (2)  No.  1,  and  the  second  line  gives 

(9). 


, dV'  , , . 

0 = M being  evanescent) 
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From  the  nature  of  the  question  de  mmimo  just  resolved,  there  can 
be  little  doubt  but  that  the  equations  (2')  and  (9)  will  suffice  for  the 
complete  determination  of  V,  where  V"  and  are  both  given.  Eut 
as  the  truth  of  this  will  be  of  consequence  in  what  follows,  we  will, 
before  proceeding  farther,  give  a demonstration  of  it;  and  the  more 
willingly  because  it  is  simple  and  very  general. 

4.  Now  since  in  the  expression  (5)  ti  is  always  positive,  every  one 
of  the  elements  of  this  expression  will  therefore  be  positive ; and  as 
moreover  V"  and  Fj'  are  given,  there  must  necessarily  exist  a function 
Fo  which  will  render  the  quantity  (5)  a proper  minimum.  But  it 
follows,  from  the  principles  of  the  Calculus  of  Variations,  that  this 
function  F’o,  whatever  it  may  be,  must  moreover  satisfy  the  equations 
(2')  and  (9).  If  then  there  exists  any  other  function  F,  which  satisfies 
the  last-named  equations,  and  the  given  values  of  V"  and  F,',  it  is  easy 
to  perceive  that  the  function 

F=  K + A{V,  - Fo) 

will  do  so  likewise,  whatever  the  value  of  the  arbitrary  constant  quan- 
tity A may  be.  Suppose  therefore  that  A originally  equal  to  zero 
is  augmented  successively  by  the  infinitely  small  increments  ^A,  then 
the  corresponding  increment  of  V will  be 

and  the  quantity  (5)  will  remain  constantly  equal  to  its  minimum 
value,  however  great  A may  become,  seeing  that  by  what  precedes 
the  variation  of  this  quantity  must  be  equal  to  zero  whatever  the 
variation  of  F may  be,  provided  the  foregoing  conditions  are  all  satis- 
fied. If  then,  besides  F^,  there  exists  another  function  satisfying 
them  all,  we  might  give  to  the  partial  differentials  of  F,  any  values 
however  great,  by  augmenting  the  quantity  A sufficiently,  and  thus 
cause  the  quantity  (5)  to  exceed  any  finite  positive  one,  contrary  to 
what  has  just  been  proved.  Hence  no  such  value  as  F^  exists. 

We  thus  see  that  when  F"  and  F^  are  both  given,  there  is  one 
and  only  one  way  of  satisfying  simultaneously  the  partial  differential 
equation  (2),  and  the  condition  (9). 
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5.  Again,  it  is  clear  that  the  condition  (4)  is  satisfied  for  the  whole 
of  F'a';  and  it  has  before  been  observed  (No.  2.)  that  when  V is  deter- 
mined by  the  formula  (1),  it  may  always  be  expanded  in  a series  of 
the  form 

V'  = A + Cu^  + &c. 

Hence  the  right  side  of  the  equation  (9)  is  a quantity  of  the  order 
; and  u'  being  evanescent,  this  equation  will  then  evidently  be 
satisfied,  provided  we  suppose,  as  we  shall  in  what  follows,  that 

« - 5 -f  1 is  positive. 

If  now  we  could  by  any  means  determine  the  values  of  V and 
belonging  to  the  expression  (1),  the  value  of  F'  would  be  had 
without  integration  by  simply  satisfying  (2')  and  (9),  as  is  evident  from 
what  precedes.  But  by  supposing  all  the  constant  quantities  «i,  «2j  «3 
a,  and  h infinite,  it  is  clear  that  we  shall  have 

0 = F", 

and  then  we  have  only  to  find  F"/,  and  thence  deduce  the  general 
value  of  F. 


6.  For  this  purpose  let  us  consider  the  quantity 


fdxidx2...dxsduu”  " 


dFdU 

dxi  dxi 


dFdU\ 

dx2  dx2  dxs  dxs  du  du  j ’ 


(10) 


the  limits  of  the  multiple  integral  being  the  same  as  those  of  the 
expression  (5),  and  U being  a function  of  x^,  x^, Xs  and  ?i,  satis- 
fying the  condition  0 = U"  when  a^, a,  and  h are  infinite. 


But  the  method  of  integration  by  parts  reduces  the  quantity  (10)  to 


— fdx^  dx2 dx, 


dU' 

du 


u 


. F' 


( d^  U 

— fdx^dx2.-‘dXsduu'~^F^i'^^  + 


d^U  ^ d^U  ^ n — s d U\ 
dti^  u du) 


.(11) 


since  0 = F'’\  and  as  we  have  likewise  0 = U",  the  same  quantity  (10) 
may  also  be  put  under  the  form 
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dV 

jdxxdXi dxg-^  . U' 


— jdxi  dxi...dxsdu  u’'~’ . U |s 


.,d^V  d^V  n — sdV\ 

5+1 • I I 


dxr  ^ dir 


u 


du  j’ 


.(12). 


Supposing  therefore  that  U like  V also  satisfies  the  equation  (2'), 
each  of  the  expressions  (11)  and  (12)  will  be  reduced  to  its  upper  line, 
and  we  shall  get  by  equating  these  two  forms  of  the  same  quantity : 


Idxi  dx2...dXi 


dU' 

du 


V'  = jdxi  dXi...dx, 


dV 

du 


u 


'n-s 


U'\ 


the  quantities  bearing  an  accent  belonging,  as  was  before  explained,  to 
one  of  the  extreme  limits. 


Because  V satisfies  the  condition  (9),  the  equation  immediately  pre- 
ceding may  be  written 


idx.  dx2...dxs^~^u 
du 


dV 

' " V = f d XidXi... dxs  u' " " ’ LV- 


If  now  we  give  to  the  general  function  U the  particular  value 

u = {(a’i  - x/y  + {x.,  - X2y  4- + (^,  - Xs'y  + u^  ; 

which  is  admissible,  since  it  satisfies  for  V to  the  equation  (2),  and  gives 
U”  = 0,  the  last  formula  will  become 


dxi  dx-i 


m'""* 


dF: 

du 


n^\ 

— x/y  + {x2  - X2y  + + {xs  - xjy  + 


=/ 


dxi  dXi dx, . (1  — w)  m'"  V' 


+ (^2  - ^2")*  + + (^*  - x',y  + 


n-H 

n\  2 


.(13); 


in  which  expression  ii  must  be  regarded  as  an  evanescent  positive 
quantity. 


In  order  now  to  effect  the  integrations  indicated  in  the  second 
member  of  this  equation,  let  us  make 


r 


ATTRACTIONS  OF  ELLIPSOIDS  OF  VARIABLE  DENSITIES.  405 


X,  - x"  = n'p  cos  9 1 ; x-i—x"  = u'p  sin  9i  cos  02 ; x^—Xz—u'p  sin  0i  sin  9^  cos  9^,  kc. 
until  we  arrive  at  the  two  last,  viz,, 


- x'^^i  = u'p  sin  9^  sin  9^ sin  9s-2  cos  0,_i, 

X,  — x"  = w'jo  sin  0^  sin  02 sin  0,_2  sin  0i_i ; 


?/'  being,  as  before,  a vanishing  quantity. 


Then  by  the  ordinary  formulee  for  the  transformation  of  multiple 
integrals  we  get 

dx^  dx^..., dxs  — i^'d'^  sin  01®"^  sin  02®'® sin0^_gC?/)  </0i  d9^...d9s-\, 

and  the  second  number  of  the  equation  (13)  by  substitution  will  become 

fdp  d9x  d9i d9,_xd~'  sin  0i®'®  sin  02®“® sin  0,_2 . (1  - n)F" 

{i+pi 


.(14). 


But  since  u'  is  evanescent,  we  shall  have  p infinite,  whenever  Xi,  Xi,...x, 
differ  sensibly  from  x",  x”,...x”\  and  as  moreover  w — # + 1 is  positive, 
it  is  easy  to  perceive  that  we  may  neglect  all  the  parts  of  the  last 
integral  for  which  these  differences  are  sensible.  Hence  V may  be 
replaced  with  the  constant  value  VI  in  which  we  have  generally 


Again,  because  the  integrals  in  (14)  ought  to  be  taken  from  0,_i  = O 
to  0,_i  = 27r,  and  afterwards  from  0^  = 0 to  0^=7r,  whatever  whole  number 
less  than  5—1  may  be  represented  by  r,  we  easily  obtain  by  means  of 
the  well  known  function  Gamma: 


/ sin  01®"®  sin  02®"®  sin  03®"^ sin  0s_2 d9^d92...d9s-x 

and  as  by  the  aid  of  the  same  function  we  readily  get 


3 g2 
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the  integral  (14)  will  in  consequence  become 


— 2 TT^ . r 


w - « + 1 
~2 


( 


n — I 


V 

^ 0 J 


and  thus  the  equation  (13)  will  take  the  form 

dV 

dx,  dXi dxsU’'~‘  i 

r du 

J \{x,-x; 


— 2 . r 


w — ^ + 1 


J + {x,  - xlj  +....  + (^,  _ xlj  + 


VL 


In  this  equation  V is  supposed  to  be  such  a function  of  x^,  x.^ x, 

and  u,  that  the  equation  (2)  and  condition  (9)  are  both  satisfied.  More- 
over V = 0,  and  Fo'  is  the  particular  value  of  V for  which 

x^  = x”’,  Xi  = x”  \ Xs  = x",  and  u = 0. 


Let  us  now  make,  for  abridgment, 


P = 


ST’  “ = 


(A), 


and  afterwards  change  x into  x',  and  x"  into  x in  the  expression  im- 
mediately preceding,  there  will  then  result 


/ 


dxi  dx-l dxlP( 


— 27r^  r 


n—s-\-\ 

2 


{ {xl-x^  + {xi-x.^^  + . . . + (a:/ -a:,)®  + u^\  ^ 


.r'...(i5), 


P being  what  P becomes  by  changing  generally  x^  into  ar/,  the  unit 
attached  to  the  foot  of  P'  indicating,  as  before,  that  the  multiple 
integral  comprises  only  the  values  admitted  by  the  condition  (a),  and 
V being  what  V becomes  when  we  make  u = 0. 

The  equation  just  given  supposes  ti  evanescent;  but  if  we  were  to 
replace  u with  the  general  value  in  the  first  member,  and  make  a 
corresponding  change  in  the  second  by  replacing  V with  the  general 
value  F',  this  equation  would  still  be  correct,  and  we  should  thus  have 
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/ 


dx^dxi dxs  P^' 

{ {x,' - + (x/ -Xif+...  H-  (x/ - Xsf +u^\~ 


r...(l6). 


For  under  the  present  form  both  its  members  evidently  satisfy  the 
equation  (2),  the  condition  (9),  and  give  V = 0.  Moreover,  when  the 
condition  (3)  is  satisfied,  the  same  members  are  equal  in  consequence 
of  (15).  Hence  by  what  has  before  been  proved  (No.  4),  they  are 
necessarily  equal  in  general. 

By  comparing  the  equation  (16)  with  the  formula  (1),  it  will  become 
evident,  that  whenever  we  can  by  any  means  obtain  a value  of  V satis- 
fying the  foregoing  conditions,  we  shall  always  be  able  to  assign  a value 
of  p which  substituted  in  (1)  shall  reproduce  this  value  of  V.  In  fact, 
by  omitting  the  unit  at  the  foot  of  P,  which  only  serves  to  indicate 

the  limits  of  the  integral,  we  readily  see  that  the  required  value  of  p is 


7.  The  foregoing  results  being  obtained,  it  will  now  be  convenient 
to  introduce  other  independent  variables  in  the  place  of  the  original 
ones,  such  that 

^1  — ^2  ^2^25... ^ ^ “ /iVf 

«i,  being  functions  of  /i,  one  of  the  new  independent  variables, 

determined  by 

di  = + /^^ 

and  V a function  of  the  remaining  new  variables,  ^i,  ^2,  iz, satis- 

fying the  equation 

1 = + + + ?/; 

a/,  «/,  d^, a/  being  the  same  constant  quantities  as  in  the  equation 

(a),  No  1.  Moreover,  «i,  will  take  the  values  belonging  to 

the  extreme  limit  before  marked  with  two  accents,  by  simply  assigning 
to  /^  an  infinite  value. 


408 


Mr  green,  on  THE  DETERMINATION  OF  THE 


The  easiest  way  of  transforming  the  equation  (2)  will  be  to  remark, 
that  it  is  the  general  one  which  presents  itself  when  we  apply  the 
Calculus  of  V ariations  to  the  quantity  (5),  in  order  to  render  it  ' a 
minimum.  We  have  therefore  in  the  first  place 

and  by  the  ordinary  formula  for  the  transformation  of  multiple  integrals, 
dx,dx, dx,du  = d^,dl,...dldh. 

But  since  1 — , 

Ur  «r 

the  expression  (5)  after  substitution  will  become 

/f/^if/^2 d^,dh  a,  «3 


Applying  now  the  method  of  integration  by  parts  to  the  variation  of 
this  quantity,  by  reduction,  we  get  for  the  equivalent  of  (2)  the  equation 

d^r  ( dV 


0 = -77^  + 


dJr 


(--a 


HI  * ('  -ifl  s ™ . I— s|l  - 


^ a 

+ X 2-^-^  - 

Ur  Uyd^f 


cCd^^ 

Ir  1.  d^V 

ar  crr>  d^rd^r' 


di: 

(2") 


+ Ul. 


Uyd^r 


Clr 


afd^r  ’ 


where  the  finite  integrals  are  all  supposed  taken  from  r = 1 to  ;’  = .?  + !, 
and  from  r'  = 1 to  r'  = s + 1. 


The  last  equation  may  be  put  under  the  abridged  form, 

d^r  ( 

provided  we  have  generally 


^ ~ dk  I”  ^ a\}  hdh  ^ ^ 
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J 2 


coefficient  of 


in  vr=  ^ ^ I?  + ^ 


«r 


Cly- 


coefficient  of  -/>.  in  v ^ > 


d^,d^, 

dV 


a?  a/ 


coefficient  of  in  v?^=-^|-w  + 2 

d/fTy  d/Y  ^ d d/y  j 


Moreover,  when  we  employ  the  new  variables 


df^  _ /,  \ h'ir  dV  dV\ 

du  ~ ^ ^ ^ . ) • Y'  ^ 2 j’ 


and  therefore  the  condition  (9)  in  like  manner  will  become 


0 = + 


1-S 


a?  j ’ \ a,?  d^r 


dkj 


(9'); 


where  the  values  of  the  variables  ^2, must  be  such  as  satisfy 

the  equation  1;“  = 0,  whatever  h may  be ; and  as  w - # + 1 is  positive,  it 
is  clear  that  this  condition  will  always  be  satisfied,  provided  the  partial 
differentials  of  F'  relative  to  the  new  variables  are  all  finite. 


8.  Let  us  now  try  whether  it  is  possible  to  satisfy  the  equation 
(2"')  by  means  of  a function  of  the  form 

Hep {ft); 

H depending  on  the  variable  h only,  and  0 being  a rational  and  entire 
function  of  ^2, 6 of  the  degree  7,  and  quite  independent  of  //. 


By  substituting  this  value  of  V in  (2"')  and  making 


^ ~ d/i^  ^ 


n 


d^\ 

-2^1 

a,') 


we  readily  get 


dH 
hdh 

0 = v0  - /f0  .••• 


+ kH. 


(17) , 

(18) ; 


where,  in  virtue  of  (17)  k must  necessarily  be  a function  of  h only; 
and  as  the  required  value  of  ep,  if  it  exist,  must  be  independent  of  //, 
we  have,  by  making  h = 0 in  the  equation  immediately  preceding, 

0 = v'0  - h<p (19) ; 

/i'o  being  the  value  k,  and  v'  0 that  of  v 0 when  /^  = 0. 
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We  shall  demonstrate  almost  immediately  that  every  function  <p  of 
the  form  (20),  No.  9,  which  satisfies  the  equation  (19),  and  which  there- 
fore is  independent  of  /i,  will  likewise  satisfy  the  equation  (18);  and 
the  corresponding  value  of  k obtained  from  the  latter  being  substituted 
in  the  ordinary  differential  equation  (17),  we  shall  only  have  to  integrate 
this  last  in  order  to  have  a proper  value  of  V'. 

9.  To  satisfy  the  equation  (19)  let  us  assume 

= &c (20); 

F being  the  characteristic  of  a rational  and  entire  function  of  the 
degree  27',  and  the  most  general  of  its  kind,  and  he,,  designating 

the  variables  in  (p  which  are  affected  with  odd  exponents  only ; so  that 
if  their  number  be  v we  shall  have 


7 = 27'  + y, 

the  remaining  variables  having  none  but  even  exponents.  Then  it  is 
easy  to  perceive,  that  after  substitution  the  second  member  of  the 
equation  (19)  will  be  precisely  of  the  same  form  as  the  assumed  value 
of  (p,  and  by  equating  separately  to  zero  the  coefficients  of  the  various 

powers  and  products  of  we  shall  obtain  just  the  same 

number  of  linear  algebraic  equations  as  there  are  coefficients  in  0,  and 
consequently  be  enabled  to  determine  the  ratios  of  these  coefficients 
together  with  the  constant  quantity  h. 


In  fact,  by  writing  the  foregoing  value  of  (p  under  the  form 

= ..f,"'?.-* f." (ao'); 

and  proceeding  as  above  described,  the  coefficient  of  ^i"i 
will  give  the  general  equation 

\J  = \^\  ^ -T  /^OJ  mg 


I’ 


-C>  V'  ^ ) /f  (01\ 


{mr+2)  {mr+  1) 


tm.>  nigj • • •m^+2> 


ATTRACTIONS  OF  ELLIPSOIDS  OF  VARIABLE  DENSITIES. 


411 


the  double  finite  integral  comprising  all  the  values  of  r and  r , except 
those  in  which  r = r , and  consequently  containing  when  completely 
expanded  s {s  — 1)  terms. 

For  the  terms  of  the  highest  degree  7 and  of  which  the  number  is 

7 ' + 1 . 7'  + 2 7'  + ^ — 1 _ 

1.2. 3 1 ~ 

the  last  line  of  the  expression  (21)  evidently  vanishes,  and  thus  we 
obtain  N distinct  linear  equations  between  the  coefficients  of  the  degree 
7 in  <p  and  kg. 

Moreover,  from  the  form  of  these  equations  it  is  evident  that  we 
may  obtain  by  elimination  one  equation  in  kg  of  the  degree  JV,  of 
which  each  of  the  roots  will  give  a distinct  value  of  the  function 
(p^'y\  having  one  arbitrary  constant  for  factor;  the  homogeneous  function 
being  composed  of  all  the  terms  of  the  highest  degree,  7 in  cp. 
But  the  coefficients  of  and  kg  being  known,  we  may  thence  easily 
deduce  all  the  remaining  coefficients  in  0,  by  means  of  the  formula  (21). 

Now,  since  the  N'  linear  equations  have  no  terms  except  those  of 
which  the  coefficients  of  0^^^  are  factors,  it  follows  that  if  kg  were  taken 
at  will,  the  resulting  values  of  all  these  coefficients  would  be  equal  to 
zero.  If  however  we  obtain  the  values  of  iV  — 1 of  the  coefficients 
in  terms  of  the  remaining  one  A from  iV-  1 of  the  equations,  by  the 
ordinary  formulas,  and  substitute  these  in  the  remaining  equation,  Ave 
shall  get  a result  of  the  form 

K.A  = 0, 

where  is  a function  of  kg  of  the  degree  JV.  We  shall  thus  have 
only  two  cases  to  consider:  First,  that  in  which  A = 0,  and  consequently 
also  all  the  other  coefficients  of  0^^^  together  with  the  remaining  ones 
in  0,  as  will  be  evident  from  the  formulae  (21).  Hence,  in  this  case 

0 = 0: 

Secondly,  that  in  Avhich  kg  is  one  of  the  JV  roots  of  0 = IC,  as  for 
instance,  kg  in  this  case  all  the  coefficients  of  0 will  become  multiples 
of  A,  and  we  shall  have 

VoL.  V.  Part  III.  sH 
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cp  = A (pi'. 

(pi  being  a determinate  function  of  ^i,  ^2, 

We  thus  see  that  when  we  consider  functions  of  the  form  (20) 
only,  the  most  general  solution  that  the  equation 

0 = S7'(p  - ho  (p (19') 

admits  is  

or,  <p  = 0;  or,  (p  — a(P', 

a being  a quantity  independent  of  ^1,  ^2, ^s,  and  (p  any  function 

which  satisfies  for  (p  to  the  equation  (19').  But  by  affecting  both  sides 
of  the  equation 

0 = v'  0 - ^0'  0 

with  the  symbol  v,  we  get 

0 = V . v'  0 - h„' . V 0 ; 

and  we  shall  afterwards  prove  the  operations  indicated  by  y and  v' 
to  be  such,  that  whatever  0 may  be, 

V v'  0 = v'  V 0. 

Hence,  the  last  equation  becomes 

v'  ( V 0)  - ^0'  V 0 ; 

and  as  V 0 like  0 is  of  the  form  (20),  it  follows  from  what  has  just 
been  shewn,  that 

either  0 = v 0,  or,  v 0 = a 0, 

a being  a quantity  independent  of  ^i,  ^2, 

The  first  is  inadmissible,  since  it  would  give  (p  - O',  therefore  when 
(p  satisfies  (19'),  we  have 

^ (p'  = a(p,  i.  e.  0 = V 0 — «0- 

But  since  a is  independent  of  ^1,  ^2, 6,  this  last  equation  is 

evidently  identical  with  (18),  since  the  equation  (18)  merely  requires  that 
/c  should  be  independent  of  ^1,  ^2, 
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Having  thus  proved  that  every  function  of  the  form  (20)  which 
satisfies  (19)  will  likewise  satisfy  (18),  it  will  be  more  simple  to  deter- 
mine the  remaining  coefficients  of  0 from  those  of  by  means  of 
the  last  equation,  than  to  employ  the  formula  (21)  for  that  purpose. 


h 

Making  therefore  h infinite  in  (18),  and  writing  in  the  place 
of  K,  we  get 

s+i  r 


0 = 2,'+  ■ . (1  - f,=)  0-2  (22)  ?.f,- 


where  (S2)  comprises  the  — — ~ combinations  which  can  be  formed  of 
the  s indices  taken  in  pairs. 


If  now  we  substitute  the  value  of  0 before  given  (20'),  and  recol- 
lect that  for  the  terms  of  the  highest  degree  we  have  ^mr  = y,  we  shall 
readily  get 

0=(7-2w,.)  (7  + Sw?,+w-l)^,„^,  „^,...„^+(w,-f-2)  {nir-k-\)A,n^,  m^+2,...  ™^...(22), 

from  which  all  the  remaining  coefficients  in  (p  will  readily  be  deduced, 
when  those  of  the  part  0*'^'  are  known. 


10.  It  now  remains,  as  was  before  observed,  to  integrate  the  ordi- 
nary differential  equation  (17)  No.  8.  But,  by  the  known  theory  of 
linear  equations,  the  integration  of  (17)  will  always  become  more  simple 
when  we  have  a particular  value  satisfying  it,  and  fortunately  in  the 
present  case  such  a value  may  always  be  obtained  from  <p  by  simply 

changing  into  -ax  • ^^^t  if  we  represent  the  value  thus  ob- 

tained  by  Ho  we  shall  have 


dHo  ^ h 

dh  * «r\/(2«)^)’ 


and  by  a second  differentiation 


d'Ho  _ d(p  U gz-v-gs  ^'0 

dJf  d^r ' cir's/^d^)  d^r  dr  • ' d^A^r'  ’ 


3h2 
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^ ^ J 

(22)  as  before  comprising  all  the  — ~ combinations  of  the  .v  indices 

1 » ^ 

taken  in  pairs. 

Hence,  the  quantity  on  the  right  side  of  the  equation  (17),  when 
we  make  H = becomes 


^d± 


a. 


(p  Id 


d^r  \/(2a7)  d^r"  ' tt,'  2 d. 


+ K(f) 


+ 2(22) 


d'(j) 


dj^rd^/  ' a,a/'S.dJ‘  ' V " d^r  ' (ir\/{^nr‘) 

But  if  we  recollect  that  we  have  generally 


(23). 


_ a, 

■ ^/(20 


(24), 


it  is  easy  to  perceive  that  in  consequence  of  the  equation  (18)  the 
quantity  (23)  will  vanish,  and  therefore  the  foregoing  value  of  Ho 
will  always  satisfy  the  equation  (17). 


Having  thus  a particular  value  of  H,  we  immediately  get  the 
general  one  by  assuming 


H=  Hojxdh. 


In  fact,  there  thence  results 


H = KHof 


h‘~”dh 


Ho-  f/2,  «3 a,  ’ 


the  two  arbitrary  constants  which  the  general  integral  ought  to  con- 
tain being  K,  and  that  which  enters  implicitly  into  the  indefinite  in- 
tegral. But  the  condition  0 = V requires  that  H should  vanish  when 
h is  infinite,  and  consequently  the  particular  value  adapted  to  the 
present  investigation  is 


TJ  TT  r 

■Ho  — /v  . Ho  / rj-9 

Hdux,  a.,... 


a. 


11.  The  values  of  0 and  H being  known,  we  may  readily  find 
the  corresponding  values  of  V and  p.  For  we  have  immediately 
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V=H<p  = K<pH,  f 

CO 


/^  ' ’‘d/i 

HoUl,  (Is 


(26), 


and  as  the  function  (p  is  rational  and  entire,  and  the  partial  differen- 
tial of  V relative  to  /^  is  finite,  it  follows  that  all  the  partial  differ- 
entials of  V are  finite;  and  consequently,  by  what  precedes  (No.  7.) 
the  condition  (9')  is  satisfied  by  the  foregoing  value  of  V,  as  Avell  as 
the  equation  (2)  and  condition  0 = V".  Hence  the  equations  {h)  and 
(c)  No.  6 will  give,  since 


dV 

du 


V 


a?  ' d^r 


dV\ 
dh  1 ’ 


and  /^  must  be  supposed  equal  to  zero  in  these  equations 


P 


“2 


dV 

dh 


(where  h = 0) : 


since  where  /^  = 0,  «,  = «/ ; and  therefore. 


1 - = 1 - 


ar 


If  now  we  substitute  for  F'  its  value  (26),  and  recollect  that  — a ■ + 1 is 

always  positive,  we  get 


- r 


u — 1 


P = 


n—s  + 1 


0 


K 


«/,  r//. 


(27), 


since  it  is  clear  from  the  form  of  that  this  quantity  may  always 
be  expanded  in  a series  of  the  entire  powers  of  In  the  preceding 
expression,  (27),  indicates  the  value  of  when  h = 0,  and  (p' 

the  corresponding  value  of  0 or  that  which  would  be  obtained  by 

simply  changing  the  unaccented  letter  into  the  accented 

ones  ^i',  Is  deduced  from 

('y)  ~ ^1  ? ^2  ^^2  ?2  ? ~ • 
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It  will  now  be  easy  to  obtain  the  value  of  V corresponding  to 


/ 

P = 


§5)  < 


(28), 


without  integrating  the  formula  (1)  No  1,  where  F is  the  character- 
istic of  any  rational  and  entire  function.  In  fact  it  is  easy  to  see  from 
what  precedes  (No.  9),  that  we  may  always  expand  in  a finite  series 
of  the  form 


F{x(t  xi xl)  = &o^o'  + ^10/  + ^2  0/  + ^3  03'  + &C. 

after  A’/,  x-l,  ho.,  have  been  replaced  with  their  values  (7).  Hence,  we 
immediately  get 

= y"-*-'  . {^o0n'  + ^10/  + ^202'  + &C.} (29). 

By  comparing  the  formulas  (26)  and  (27)  it  is  clear  that  any  term, 
as  hr<^y  for  instance,  of  the  series  entering  into  p\  will  have  for  cor- 
responding term  in  the  required  value  of  V,  the  quantity 


27r*r 


« — 5 -1- 1 


W — 1' 


«/  «/ a/  . ir0r  J 


h^-”dh 


«i  tta a. 


.(30): 


Ho  being  a particular  value  of  H satisfying  the  equation  (17),  and 
immediately  deducible  from  0 by  the  method  before  explained. 


12.  All  that  now  remains,  is  to  demonstrate  that 

V'v0  = vv'0 (31), 

whatever  0 may  be.  For  this  purpose  let  us  here  resume  the  value 
of  A0,  as  immediately  deduced  from  the  equation  (2")  No.  7,  viz. 

a^a’^ 

-H  X 

a*«J  a,  «'<*{ 


• (32), 
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where  for  simplicity  the  indices  at  the  foot  of  the  letters  ^ and  a have 
been  omitted,  and  their  accents  transferred  to  the  letters  themselves. 
Moreover  all  the  finite  integi'als  are  supposed  taken  from  1 to  ^+1. 

By  making  /^  = 0 in  the  last  expression  we  immediately  get  v'<p, 
and  if  for  a moment,  to  prevent  ambiguity,  we  write  hr  in  the  place 
of  the  original  and  omit  the  lower  indices  as  before,  we  obtain 

V>=(l-2r)2j.^+  - «-l)2|^ (3S); 

wheTe  to  avoid  all  risk  of  confusion  r has  been  changed  into  r”,  and 
the  double  accent  of  this  index  transferred  to  the  letters  ^ and  h 
themselves. 


We  will  now  conceive  the  expression  (32)  to  be  written  in  the 
abridged  form 

= Vl0  + + /^W40  - 

the  order  of  the  terms  remaining  unchanged. 

If  then  we  recollect  that  the  accents  have  no  other  office  to  per- 
form than  to  keep  the  various  finite  integrations  quite  distinct,  and 
consequently  that  in  the  final  results  they  may  be  permuted  in  any 
way  at  will,  we  shall  readily  get 


(1-2^)1422: 


i 1 

1 

\d^¥ 

V Vi0  — ViV  0 — 

+22^x2^  -22lx2-f/^,. 

(t  (2)  ¥ 


(1) 


t (3) 


+ 2(^-w-1)|2:£x  s 


-2^x2-^  1 

V V20  — V2V  0 = 

(6)  "^^'(7)1 


+22-^x2:|x 


rt’  a^b"^' 


(8) 


tti  ¥d^ 


(9) 
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+ 2(^ 


V3V'0  - V'V30  = 


(1  - 


V V4^  V4V  0 — 


VsV  0 — V V50  — 


all  the  finite  integrals  being  taken  from  r = 1 to  ;•  = a-  + 1,  and  from 
r'  = 1 to  r'  = 5 + 1. 

In  order  to  obtain  tbe  required  value 


v'v0  - vv'0, 


it  is  clear  that  we  shall  only  have  to  add  the  first  of  the  five  pi’eceding 
quantities  to  the  sum  of  the  four  following  ones  multiplied  by  //*,  and 
to  render  this  more  easy,  we  have  appended  to  each  of  the  terms  in  the 
])receding  quantities  a number  inclosed  in  a small  parentbesis. 

Now  since  the  accents  may  be  permuted  at  will,  and  we  have  like- 
wise cr  = + h\  it  is  easy  to  see  that  the  terms  marked  (1),  (6)  and 

(12)  mutually  destroy  each  other.  In  like  manner,  (2),  (3),  (7)  and 

(18)  mutually  destroy  each  other;  the  same  may  evidently  be  said  of 

(13)  and  (1 6),  of  (15)  and  (17),  of  (9)  and  (19),  and  of  (8)  and  (14). 

Moreover  the  four  quantities  (4),  (5),  (10)  and  (11)  will  do  so  likewi.se, 
and  consequently,  we  have 


v'v0  - V v'0  = 0. 

Hence  the  truth  of  the  equation  (31)  is  manifest. 


ATTRACTIONS  OF  ELLIPSOIDS  OF  VARIABLE  DENSITIES. 


419 


Application  of  the  preceding  General  Theory  to  the  Determination  of  the 

Attractions  of  Ellipsoids. 

13.  Suppose  it  is  required  to  determine  the  attractions  exerted  by 
an  ellipsoid  whose  semi-axes  are  a',  b',  c'  whether  the  attracted  point 
j)  is  situated  within  the  ellipsoid  or  not,  the  law  of  the  attraction  being 
inversely  as  the  power  of  the  distance.  Then  it  is  well  known 
that  the  required  attractions  may  always  be  deduced  from  the  function 

p dx'  dy  d% 

{{x~xy  + {y-  y'f  f 

p being  the  density  of  the  element  dx'  dy'  d%'  of  the  ellipsoid,  and 
X,  y,  z being  the  rectangular  co-ordinates  of  p. 

We  may  avoid  the  breach  of  the  law  of  continuity  which  takes 
place  in  the  value  of  V,  when  the  point  p passes  from  the  interior  of 
the  ellipsoid  into  the  exterior  space,  by  adding  the  positive  quantity 
to  that  inclosed  in  the  braces,  and  may  afterwards  suppose  ti  eva- 
nescent in  the  final  result.  Let  us  therefore  now  consider  the  function, 

p dx'  dy  dz' 

{ (^  - x’y  + (y-  y'f  + (z  — z'f  + u^]  * 

this  triple  integral  like  the  preceding  including  all  the  values  of  x',  y',  z', 
admitted  by  the  condition 

X'^  y'2  g,/2 

— -F  — + — Z 1 
o',  -t-  ^ 


^=/ 


If  now  we  suppose  the  density  p’  is  of  the  form 


'2s  n'-2 


S''. *') 


(84), 


which  will  simplify  f {x',  y,  z')  when  p is  constant  and  n'  = 2,  and  then 
compare  this  value  with  the  one  immediately  deducible  from  the  general 
expression  (28)  by  supposing  for  a moment  »'  = w,  viz. 


p — ^^2  ^'2  j T {x , y , z). 


VoL,  V.  Part  III. 
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we  see  that  the  function  f will  always  be  two  degrees  higher  tlmn 
But  since  our  formula  become  more  complicated  in  proportion  as  the 
degree  of  F is  higher,  it  will  be  simpler  to  determine  the  differentials 
of  V,  because  for  these  differentials  the  degree  of  F and  f is  the  same. 
I..et  us  therefore  make 


A = 


1 

1 - n' 


dV 

dx 


/- 


p {x  — x)  dx  dy  d! 


n'  + 1 > 


then  this  quantity  naturally  divides  itself  into  two  parts,  such  that 

A = xA  + A\ 


where  A — + J- 
and  A"=  — J 


p'  dx'  dy  dz 


n'  + I ’ 


{{x-x'f-{-{y-y'f  + {z-zf  + u^}  « 
x'  p dx'  dy  dz 

{ {x  - x'f  + {y-  y'f  A{z-  z’f  + \ 


n*+  I * 


By  comparing  these  with  the  general  formula  (1),  it  is  clear  that 
u — 1 = n'  + 1,  and  consequently  n = n + 2.  In  this  way  the  expression 
(28)  gives 


r'ii  2.'2\  ”'-2 


which  coincides  with  (34)  by  supposing  F=f. 


'fhe  simplest  case  of  the  present  theory  is  where  f{x!,  y\z)  = \,  and 
then  by  No  11,  we  have  0o'=  1 and  Z»o  = 1,  when  A is  the  quantity 
required,  and  as  the  general  series  (29),  No  11,  then  reduces  itself  to 
its  first  term,  we  immediately  obtain  from  the  formula  (30),  the  value 
of  A following. 


A = 


M r (I 
V2 


n + 1 
2 


ahc 


(35), 


because  in  the  present  case  =1,  s = 3,  and  n = n -|-  2. 


Again,  the  same  general  theory  being  applied  to  the  value  of  A" 
given  above,  we  get 
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F {x,  y,  z')  = - x'f(x,  y',  z')  =-x'  (when  /=  1), 

and  hence  by  Noll,  F{x,y,z')  = —a^.  In  this  way  the  series  (29) 
again  reduces  itself  to  a single  term,  in  which 


(po  = r,  and  ^0  = — (i, 

and  the  particular  value  corresponding  thereto,  by  omitting  the  super- 

1 

fluous  constant  ^ ^ ^.'2^  10)’ 

Ho  = a. 


These  substituted  in  the  general  formula  (30)  as  before,  immediately  give 


a^b'c'  j 
00 


a^bc  ’ 


and  consequently  by  reduction  since  a '^  = x. 


A=xA'  + A"  = — 


a'b'c'x  J 
00 


a^bc 


The  value  of  A just  given  belongs  to  the  density 
^ V b'^  cV 


Hence  we  immediately  obtain  without  calculation  the  corresponding 
values 

""•''S') ...... 


B = 


C = 


1 dV 
1—n'  dy 


1 dV 


n>  > rJi  dh 
— — — abcy  / — rr-  > 
+ 1^  ab^c 

(2) 

7 — — - abcz  — . 
' + 1\  abd’ 


\—n'  dz 


27r3r 


-) 


3 I 2 
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If  now  we  suppose  moreover 

\ dV  r p dx  dy  d% 


D = 


\ - n'  du 


= uf- 


n'+  1 ? 


K^-^T  + (y-yT  + («-~T+«'l  " 

the  method  before  explained  (No  11),  will  immediately  give 


n = - 


; — a b cu  / i — 

n + 1\  ahc 

I 


and  therefore  if  for  abridgment  we  make 


TT 


f r 


31  = (n'—l)  ■ 


I') 


(t) 


f j f / 

ahc, 


the  total  differential  of  V may  be  written 

,/r=M[^,:dx  + S,jd,j  + 

which  being  integrated  in  the  usual  way  by  first  supposing  h constant, 
and  then  completing  the  integral  with  a function  of  h,  to  be  after- 
wards determined  by  making  every  thing  in  V variable,  we  get 


k being  a quantity  absolutely  constant,  which  is  equal  to  zero  when 
;^'>l.  What  has  just  been  advanced  will  be  quite  clear  if  we  recollect 
that  h may  be  regarded  as  a function  of  x,  y,  z and  u,  determined  by 
the  equation 


1 ^ y 

1 = ■ /'o  - -JO  + “ 


+ + i = r + 


«'*  + h"  h'^  + c'^  + h 

seeing  that  = a'^  + 4^,  b^  = //*  + h^,  and  c*  = c~  + h\ 


(37); 


After  what  precedes,  it  seems  needless  to  enter  into  an  examination 
of  the  values  of  V'  belonging  to  other  values  of  the  density  p,  since 
it  must  be  clear  that  the  general  method  is  equally  applicable  when 
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.'=(1 

where  f is  the  characteristic  of  any  rational  and  entire  function. 


The  quantity  A before  determined  when  we  make  u = 0,  serves  to 
express  the  attraction  in  the  direction  of  the  co-ordinate  x of  an  ellipsoid 
on  any  point  p,  situated  at  will  either  within  or  without  it.  But  by 
making  u =-0  in  (37)  we  have 


1 


x~ 


+ 


y 


a'^  + h^  ' + + ' It 


+ 


(38), 


and  it  is  thence  easy  to  perceive  that  when  j)  is  within  the  ellipsoid, 
h must  constantly  remain  equal  to  zero,  and  the  equation  (38)  will  always 

o"  ' 

be  satisfied  by  the  indeterminate  positive  quantity  ^ . When  on  the 

contrary  p is  exterior  to  it,  h can  no  longer  remain  equal  to  zero,  but 
must  be  such  a function  of  x,  y,  z,  as  will  satisfy  the  equation  (38),  of 
which  the  last  term  now  evidently  vanishes  in  consequence  of  the 
numerator  o^.  Thus  the  forms  of  the  quantities  A,  B,  C,  D and  V 
all  remain  unchanged,  and  the  discontinuity  in  each  of  them  falls  upon 
the  quantity  h. 


To  compare  the  value  of  A here  found  with  that  obtained  by  the 
ordinary  methods,  we  shall  simply  have  to  make  n'  = 2 in  the  expression 

(36),  recollecting  that  r(l)  = 1,  and  F TT.  In  this  way 

A = — ^ira'b'c'x  I = — ^iira'b'c'  X f 


dbc 


= + 4ia'b'c'x  f ~~  = ‘^iwa'b'c  f - 

Ja  dbc  Ja 


da 


a'^  + b'^)  {d-a’^  + c-) 


But  the  last  quantity  may  easily  be  put  under  the  form  of  a definite 
integral,  by  writing  ^ in  the  place  of  a under  the  sign  of  integration, 

and  again  inverting  the  limits.  Thus  there  will  result 
^Trdb'c  r'  v^dv 


A = 


/ 


— ct^ 

1 + 


a 


«" 
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which  agrees  with  the  ordinary  formula,  since  the  mass  of  the  ellipsoid 
^Trdh'c 


IS 


3 


and  «-  = + h^. 


Kxamintttion  of  a particula?'  Case  of  the  General  Theory  exposed  in  the 

former  Part  of  this  Paper. 

14.  There  is  a particular  case  of  the  general  theory  first  considered, 
which  merits  notice,  in  consequence  of  the  simplicity  of  the  results  to 
wliich  it  leads.  The  case  in  question  is  that  where  we  have  generally 
whatever  r may  be 


a;  = a. 


Then  the  equation  (19)  which  serves  to  determine  (p,  becomes  by 
supposing  ko  = k . a'^ 


.(39). 


0 = + (.9-w-l) 

If  )iow  we  employ  a transformation  similar  to  that  used  in  obtaining 
the  formula  (14),  No  6,  by  making 

^i  = P cos  9,,  ^2  = p sin  9i  cos  ^3  = p sin  sin  G,>  cos  $3,  &c. 
and  then  conceive  the  equation  (39)  deduced  from  the  condition  that 

{2,-  (^)’  + 

must  be  a minimum  (vide  No  8),  we  shall  have 

r/^,  r/^2 d^,  = .sin0,*“^  sin^z*"* sin  6, .2  dpdO^  dGi d0,_^, 

+ l 4-  - S*  

\dfj  \dp}  d ‘ sin  0,®  sin  02* sin0*,_,’ 

and  1 - 2^/  = 1 - p\ 

Proceeding  now  in  the  manner  before  explained,  (No  8),  we  obtain 
for  the  equivalent  of  (39)  by  reduction 

d^(p  , _ . _ 1 \ cos  Or  d<p 

d^(p  .V  - 1 — np^  d^  1 ^ ^ dG?  ^ ^ ' sin  0,  dGr  h 


0 = 


+ 


dp^  p{l-p^)  'dp  |0®  ‘ sin  0,®  sin  02® sin  0®^_, 


1 - 


10. ..(40). 
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But  this  equation  may  be  satisfied  by  a function  of  the  form 

<P  — jP0i0203 0s-i; 

P being  a function  of  p only,  and  afterwards  generally  0,  a function 
of  6r  only.  In  fact,  if  we  substitute  this  value  of  (p  in  (40),  and  then 
divide  the  result  by  (p,  it  is  clear  that  it  will  be  satisfied  by  the  system 

0,,_ir/0\_i  “ 

(l‘'Os-2  .J  COS  9s-2  (lOs_2  ^s-1  ^'A'l  ^ 

Os-2(1&\^2  ‘ sin  Qs-idds-2  sin  ~ ^ ’ 

^"0.-3  2 COS  </0,-3  X,_2  _ 

^,-id0'\_2  ■ sin  0s_3  0,_3G?0,_3  sine^_3  ~ 

&c.  &c.  &c.  &c. 


combined  with  the  following  equation, 

d^P  s—\  — np^  dP  \ik 

pd^^  ^ ■ Pd^  7 ~ ^ 

where  k,  Xi,  Xa,  X,,  Ike.  are  constant  quantities. 


(42), 


In  order  to  resolve  the  system  (41),  let  us  here  consider  the  general 
type  of  the  equations  therein  contained,  viz. 


d^Qs-r  / COS0^_r  </0s_,. 

dG~,_r  ^ ^ sin  0s_/ 


+ 


( 


r+l 


\sin  0^ 


Xs  — r 1 0.'i  - r • 


Now  if  we  reflect  on  the  nature  of  the  results  obtained  in  a preceding 
part  of  this  paper,  it  will  not  be  difficult  to  see  that  0,,_,.  is  of  the  form 

Os-r  = (sin0,_,)>  = 

where  y)  is  a rational  and  entire  function  of  m = cos0s_,,  and  I a whole 
number. 


By  substituting  this  value  in  the  general  type  and  making 


we  readily  obtain 


0=(1-M’) 


d^p 


Xs_,.+  i = - i{i  ~>r  r - 2) 

(2e  + r)j(A^—  \Xs-r  + i{i  + r-  1)}  p. 


(43) 
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To  satisfy  this  equation,  let  us  assume 

p = 

Then  by  substituting  in  the  above  and  equating  separately  the  coefficients 
of  the  various  powers  of  n,  we  have  in  the  first  place  from  the  highest 

==  - e{e  + r-l) (44), 

and  afterwards  generally 

. e — i — 2t.e—i—2t—l  . 

" “ 2^  + 2x  2e  + r-2t-3  ’ 

But  the  equation  (43)  may  evidently  be  made  to  coincide  with  (44),  by 
writing  for  i,  and  for  e,  since  then  both  will  be  comprised  in 

i f^<'•>  + r-2| (45). 


Hence  we  readily  get  for  the  general  solution  of  the  system  (41), 


+ 


2 X + r — 3 


- &c.] ; 


2 . 4 X {2/^'*  + r — 3j  {2^^'■^  + r - 5} 

where  /u.  = cos0,_r,  and  represents  any  positive  integer  whatever,  pro- 
vided is  never  greater  than 


Though  we  have  thus  the  solution  of  every  equation  in  the  system 
(41),  yet  that  of  the  first  may  be  obtained  under  a simpler  form  by 
writing  therein  for  X,_i  its  value  deduced  from  (45).  We  shall 

then  immediately  perceive  that  it  is  satisfied  by 


In  consequence  of  the  formula  (45),  the  equation  (42)  becomes 
^ f/’P  s—l  — np^  dP  -t-^— 2)  k ] „ 

" = IV  " Tinyf ■ 1 — ? — ^ + w -r- 


which  is  satisfied  by  making  — X,  - (/^'* -)- 2a))  (/<’> -t- 2a)  + « — 1),  and 
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P= 


2ft)  X 2i'‘*  + 2ft)  + ^ — 2 , 

2 X 2^*"'  ^-4ft)  + ?^— 3 ^ 


+ 


2ft)  . 2ft)  — 2 X 2^'''  + 2ft)  + # — 2 . 2?'*'  + 2(0  + s—  4 
2.4x2  e’*'*  + 4ft)4-w  — 3.2*-f-4ft)  + w — 5 


.,20.-4 


— &C.} 


where  to  represents  any  whole  positive  number. 


Having  thus  determined  all  the  factors  of  tp,  it  now  only  remains 
to  deduce  the  corresponding  value  of  H.  But  Ho  the  particular  value 
satisfying  the  differential  equation  in  H,  will  be  had  from  (p  by  simply 
making  therein 

r _ CLr  a 

since  in  the  present  case  we  have  generally 


Hence,  it  is  clear  that  the  proper  values  of  6^,  62,  83,  &cc.  to  be  here 
employed  are  all  constant,  and  consequently  the  factor 

01©2  03 ©S-l 

entering  into  (p  is  likewise  constant.  Neglecting  therefore  this  factor 
as  superfluous,  we  get  for  the  particular  value  of  H, 


Ho  = Pa_; 

a' 

since  + + 


sa^  _ 
sa'^  a'^  ’ 


and  Pa  represents  what  P becomes  when  p is  changed  into  — , . 

a 


Substituting  this  value  of  Ho  in  the  equation  (25),  No  10,  there 
results  since  o*  = + P 

H=K.Pa  f , (46) 

a' 

K being  an  arbitrary  constant  quantity. 


Thus  the  complete  value  of  V for  the  particular  case  considered  in 
the  present  number  is 
VoL.  V.  Paet  III. 


sK 
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Pe.e, Q,.,.KPa_J 


k-^dJi 


■(47) 


»P„-  {d^  + df 

a' 

and  the  equation  (27),  No  11,  will  give  for  the  corresponding  value  of  p, 
n — 1'^ 


TTTITTT  ®'- 

2-«l 


where  P/,  0/,  Oj,  &c.  are  the  values  which  the  functions  P,  0i,  0a,  &c. 
take  when  we  change  the  unaccented  variables  ^i,  into  the  cor- 
responding accented  ones  and 


P.  = 


W-^-fl  .TI  — S + 3 71  — S ■{-2(0-1 

7^-^•  2e  + 2w  — 1 . w + 2i+  2w  + 1 n -\-2i  + 4-w  — 3’ 


or  the  value  of  P when  |0  = 1 ; where  as  well  as  in  wliat  follows  I 
is  written  in  the  place  of  ?***. 

The  differential  equation  which  serves  to  determine  H when  we 
introduce  a instead  of  k as  independent  variable,  may  in  the  present 
case  be  written  under  the  form 

0 = d (d  — d^)  ^ ^ + d 1 — (^  — 1) . dH 

^ ' ad  ada 

+ I?"  (z  -{- 5 — 2)  fl  ^ — (z  2 (t))  (?'  -j-  2 w + n — 1)  d^  Hf 

and  the  particular  integral  here  required  is  that  which  vanishes  when 
h is  infinite.  IMoreover  it  is  easy  to  prove,  by  expanding  in  series,  that 
this  particular  integral  is 

H = //«' A"  . 2 . 

provided  we  make  the  variable  r to  which  A"  refers,  vanish  after  all 
the  operations  have  been  effected. 


But  the  constant  //  may  be  determined  by  comparing  the  coefficient 
of  the  highest  power  of  a in  the  expansion  of  the  last  formula  with 
the  like  coefficient  in  that  of  the  expression  (46),  and  thus  we  have 


/<•'  = AV‘+2“-(_i)“ 


n + 2e  + 2w  — 1 . n + 2e  + 2£t)  + 1 w + 2/  + 4w—  3 

2.4. 6 2w 
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Hence  we  readily  get  for  the  equivalent  of  (47), 


r=pe.eo e,_,  x 


n + + 2w—l  .w  + 2i-t-2w-M w + 2^^-4w  — 3 


2 .4.6 2< 


X — 1)" a*  A*" «*’■/</« a 


s- 1 -n-2<u 


l-2r~s-2i  / 


{a^  — d^) 


In  certain  cases  the  value  of  V just  obtained  will  be  found  more 
convenient  than  the  foregoing  one  (47).  Suppose  for  instance  we  repre- 
sent the  value  of  V when  /^  = 0,  or  a = d by  V^.  Then  we  shall  hence 
get 


= P0i©2 ©s_l  X 


w -f  2?  -F  2a)  — 1 . w + 2^  + 2w  -f  1 w-F2e  + 4a)  — 3 


X ( — I)"* . A’"d^'' /‘"'dad' 


2.4.6 2a) 

S-l  ~n-2u> 


{d  — d'^) 


which  in  consequence  of  the  well  known  formula 

ra_^)r(£±f^) 


a~'^ da{d  — d'^Y^  ~ — x 

so 

by  reduction  becomes 
V,=  -PGiG^ e.-i  X 


2r 


j,  ,1  + ^ — + 4i(t}  — 1 


2 


2 


2r(a)+  i)r 


^ + 2*  -F  2a) 
2 


Kd^-” (49); 


since  in  the  formula  (^),  r ought  to  be  made  equal  to  zero  at  the  end 
of  the  process. 


By  conceiving  the  auxiliary  variable  u to  vanish,  it  will  become  clear 
from  what  has  been  advanced  in  the  preceding  number,  that  the  values 
of  the  function  V within  circular  planes  and  spheres,  are  only  particular 
cases  of  the  more  general  one,  (49),  which  answer  to  5 = 2 and  5 = 3 
respectively.  We  have  thus  by  combining  the  expressions  (48)  and 
(49),  the  means  of  determining  when  the  density  p is  given,  and 
vice  versa;  and  the  present  method  of  resolving  these  problems  seems 
more  simple  if  possible  than  that  contained  in  the  articles  (4)  and  (5) 
of  my  former  paper. 

GEORGE  GREEN. 
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